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H . ABSTRACT 

Context. Owing to the high energy required for their excitation, high-./ CO transitions are a valuable tool for the study of protostellar 
jets and outflows. However, high spectral resolution observations of high-./ CO lines, which are essential to distinguish the different 
components in the line profiles, were impossible until the start of operations of the Herschel space observatory and the Stratospheric 
Observatory for Infrared Astronomy (SOFIA). 

Aims. We present and analyze two spectrally resolved high-./ CO lines towards a protostellar outflow. We study the physical condi- 
tions, as a function of velocity, traced by such high-energy transitions in bipolar outflows. 
Methods. We selected the molecular outflow Cep E, driven by an intermediate-mass class protostar. Using the GREAT receiver on 
board SOFIA, we observed the CO (12-11) and (13-12) transitions (E„ ~ 430 and 500 K, respectively) towards one position in the 
rl^ ' blue lobe of this outflow, that had been known to display high-velocity molecular emission. 

Oh. Results. We detect the outflow emission in both transitions, up to extremely high velocities (~ 100 km s~' with respect to the systemic 

velocity). We divide the line profiles into three velocity ranges that each have interesting spectral features: standard, intermediate, and 
extremely high-velocity. One distinct bullet is detected in each of the last two. A large velocity gradient analysis for these three 
velocity ranges provides constraints on the kinetic temperature and volume density of the emitting gas, > 100 K and >10 4 cm" 3 , 
respectively. These results are in agreement with previous ISO observations and are comparable with results obtained by Herschel for 
similar objects. 

Conclusions. High-7 CO lines are a good tracer of molecular bullets in protostellar outflows. Our analysis suggests that different 
physical conditions are at work in the intermediate velocity range compared with the standard and extremely high-velocity gas. 
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1. Introduction vations of low-lying rotational transitions of the CO molecule 

-Y-! (usually up to / u = 3, corresponding to E u < 50 K). In most 

r* Observations over the past few decades have shown that during caseSi these low . energy transitions are mor e sensitive to the 

^ the star formation process, mass accretion is commonly associ- low excitation swep t- U p gas, not exclusively related with the jet 

^ . ated with mass ejection in the form of jets. The interaction be- and/or shock front Qn the other han(i higher . lying trans i tions 

. . ■ tween the jet and the parent cloud generates shock fronts, that ^ > 100 K) probe ±e wam gas unambiguous i y re i ate d with 

form large cavities called bipolar outflows. The simultaneous the jet and/or hot spots in the shocked gas Howevel - transitions 

presence of jets and bipolar outflows, as well as accretion disks wkh relatively nigh _ y rotationa i quant um numbers (e.g. 7 U > 

surrounding protostellar condensations, are commonly observed 9) fall intQ the window of the far . infrared regime that can not 

( |Arceetal. || 2007|) , making them the basic elements of the star be observed from the ground High _y co observations were pi- 



X 



formation scenario. In this picture, different excitation regimes, oneered wkh the Infmred Space observator y (IS 0), whi ch ac _ 
related to the different structures, are revealed with a diversity of ^ targeted several prot ostellar outflows (iNisini et al.l fl999. 



observational probes, such as continuum emission and molecular 2000; Gi anmnieta ij [200I . These ISO observations revealed 

transitions at different wavelengths. In the jet/outflow system, at ^ the high _ y CQ emission does indeed trace outflowing gas 

its temperature peak the shocked material is mainly radiatively aboye 1000 K> but ^ ^ resolution did not allow us 

cooled through the near-IR transitions of H 2 , whereas the low- to probe ekher the presence of differ ent spectral components, or 

lying CO transitions trace the lower density, lower temperature the physical conditions of the gas as a function its velocity In 

« 50 K) post-shock gas. In-between those two extreme regimes addilio ^ their low spatial reso i ution (,. 80 » at wavelengths cor- 

and over the 100-2000 K range, the cooling of the gas occurs responding t0 its Long Wavelength Spectrometer, LWS) compli- 

mainly through the emission of atomic and molecular lines such cates the analysis owing tQ the presence of differ ent structures 

as [OI], H 2 0, high-/ CO, OH, and H 2 transitions, which fall in witWn ±e „ beam Qnly recentlyi wkh Herschd in opera . 

the mid- to far-infrared spectral range (Nisim et aljl999|>. tion and SOF IA starting operation, we had the opportunity to 

Traditionally, the physical conditions of the gas in protostel- 
lar outflows have been studied by means of ground-based obser- 
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follow-up on these previous studies with high-resolution FIR 
spectroscopy and angular resolutions in the 15-20" range. 

In this letter, we present one of the first studies of two 
high-/ CO spectrally resolved lines in a protostellar outflow, 
as an example of the potential of SOFIA observations with 
GREAT in the context of protostellar jets/outflows. For this 
study, we selected the protostellar outflow Cep E, in the Cepheus 
molecular cloud, since it is associated with a relatively nearby 
(~730pc) interm ediate-mass class source (IRAS 2301 1 +6126: 
IChini et all 1200 lh of 100 L Q , surrounded by an envelope struc- 
ture of mass ~ 18 M . The jet/outflow syste m itself has been in- 
tensively studied a t multiple wavelengths (E isloffel et al. 1996; 
iLefloch et alJll996LlMoro-Martin et alJ200ltlSmith et al.ll2 003). 
The strong CO wings and the evidence of jet-like components 
in the spectrum of Cep E outflow make it an ideal target for the 
study of the high-/ CO emission in jet/outflow systems. 



2. Observations 

The SOFIA telescope was pointed towards the southern lobe of 
the Cep E outflow (Fig. [T}. As seen in Fig.Q] the SOFIA obser- 
vations covered the infrared knots BI and BII. The observations 
were performed during basic science flights in July 2011, with 
the Germa n Receiver for Astron omy at Terahertz Frequencies 
(GREAlQ: iHevminck et ai]|2012h . The LI band of the GREAT 
instrument was used to observe, in two separate flights, the CO 
(12-11) and (13-12) transitions at -1.4 THz and -1.5 THz, 
in the lower and upper sidebands, respectively. In Table [1] the 
line frequencies, diffraction-limited beam sizes, time ON source 
(?onX an d atmospheric opacities (r atm ) during the observations 
are presented. The front-end was connected to the AFFTS back- 
end (Klein et al. 2012), providing a total bandwidth of 1.5 GHz 
and a spectral resolution of 212 kHz, that at the observed fre- 
quencies correspond to ~ 300 km s -1 and ~0.4 km s -1 , respec- 
tively. Single pointing observations were done in double-beam 
chopped mode, with a chop throw of 60" in RA (at 1 Hz). 
The pointing was established with the optical guide cameras, 
and was stable to 4". The calibration uncertainty is within 20%. 
From observations of Jupiter and Mars, a main- beam efficiency 
of 0.5 4±0.05 in the LI band was determined (IHevminck et al.l 
1201 2l) . that - with a forward efficiency of 0.95 - was used to con- 
vert antenna to main beam temperatures, 7mb- Subsequent anal- 
ysis of the data was done with CLASS within the GILDAS soft- 
warq3, following the standard procedures of baseline subtraction 
and spectra averaging. A first-order baseline was subtracted from 
all spectra. 

Calibratio n was performed b y carefully fitting the observed 
sky emission dGuan et al.ll2012h . Unfortunately, both CO tran- 
sitions suffer from residual terrestrial atmospheric features at 
SOFIA'S flight altitude: the CO (12-11) detection bandwidth 
contains a strong water absorption (from the signal band) at 
positive velocities >60 km/s. This does not however affect the 
calibration of the pre-dominantly negative velocity emission to- 
wards the southern lobe. This is different for CO (13-12), whose 
calibration (and baselines) at extremely negative velocities (of 
lower than - 120 km s _1 ) are more uncertain owing to an image 
band atmospheric feature. However, quite comparable integrated 
intensities for both lines at small S/N values supports our confi- 
dence in the adopted calibration procedure. 



Table 1. The CO observations. 



1 GREAT is a development by the MPI fur Radioastronomie and 
the KOSMA/Universitat zu Koln, in cooperation with the MPI fur 
Sonnensystemforschung and the DLR Institut fur Planetenforschung. 

2 http://www.iram.fr/IRAMFR/GILDAS 
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<a> Tmb sca l e > at 3 km s ' spectral resolution. {b) Convolved to 20" from 
original ~ 11" resolution (IRAM-30m), Lefloch et al. in prep. 
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Fig. 1. Spitzer/IRAC band-two (4.5 fim) image of the Cep E pro- 
tostellar outflow (retrieved from the Spitzer archive). White cir- 
cles mark the diffraction-limited SOFIA beams at the frequen- 
cies of CO (12-11) and (13-12) (decreasing beam size with 
increasing / u , see Table Q]), centered at the brightest IR knot 
of the southern lobe, a =23''03 m ll s :70 6 =+61°42'06'/0 (black 
plus symbol). White letters label the two main knots inside our 
beams, BII and BI, as well as the Cep E central source (IRAS 
23011: black star). 



3. Results 

Together with the SOFIA/GREAT spectra, in Fig. [2] we also 
show a CO (2-1) line profile taken at the corresponding posi- 
tion from a Nyquist sampled map, obtained with the IRAM-30m 
telescope (Lefloch et al., in prep). The CO (2-1) spectrum was 
convolved to a beam of 20" (i.e. the SOFIA beam at the (13-12) 
transition, see Table [TJ. From the CO spectra shown in Fig. |2]at 
least three spectral features can be distinguished: two secondary 
peaks, at around -125 km s and -65 km s -1 , and a wing-like 
profile (smooth decrease of intensity towards high velocities) in 
the range from -50 km s to -20 km s . The secondary peaks 
are similar to the so-called molecular 'bullets' often observed in 
the spectrum from class protostars and likely related to the jet 
component (e.g. L1448: iBachiller et~ai]fl99 0). These 'bullets' 
were reported previously in Cep E (Lefloch et al. 1996). The 
wing-like feature is typical of the outflow phenomenon, mostly 
tracing the cavity walls crated by the jet, and usu ally referred to 
as the 'standard' high-velocity component (e.g., Bachiller et al. 
1990). On the basis of these spectral features, in Fig. |2]we di- 
vide all spectra into three velocity ranges: extremely high veloc- 
ity (EHV: -140 to -100 km s~ : ), intermediate-to-high velocity 
(IHV: -100 to -50 km s" 1 ), and standard high velocity (SHV: 
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-50 to -20 km s -1 ). In Table |2l we present the line intensities 
integrated within the defined velocity ranges. From the spectra 
in Fig. [2] one can see that the bullet-like profile in the IHV range 
is more prominent in high-/ than in low-/ CO, with the lat- 
ter still being dominated by the wing-like profile that extends 
into this velocity range. On the other hand, the SHV component 
is stronger in low-/ than in high-/ CO. From the integrated 
intensities in Table [2] we found that the (2— 1 )/( 1 2— 11) ratio is 
about 0.8, 0.4, and 2.2, for EHV, IHV, and SHV, respectively, 
indicating different excitation conditions in the different veloc- 
ity ranges. These findings underline the necessity for velocity- 
resolved spectroscopy in the excited CO transitions, to identify 
the bullets as distinct features. The significant emission of the 
high-/ CO line at high-velocity indicates that it is very likely 
related to the jet component, unlike the wing component, which 
becomes weaker in these lines. Finally, the emission at low ve- 
locities is still considerable in our high-/ CO lines. This im- 
plies that the a priori assumption made in the analysis of some 
ISO observations, namely that the excited CO em ission origi- 
nates exclusively from the high-velocity gas (e.g. Nisi ni et alj 
2000), could be incorrect. 



Table 2. Velocity ranges and integrated intensities. 




-50 
Velocity (km/s) 



Fig. 2. SOFIA/GREAT spectra of CO (12-11) and (13-12) in 
Cep E, together with the low-/ CO (2-1) transition, taken at the 
position shown in Fig.Q] The CO (13-12) and (2-1) spectra are 
taken from data with a spatial resolution of 20", while the spatial 
resolution of the CO (12-11) is 21'.'7. The two upper spectra 
have been shifted along the Y-axis. The green horizontal lines 
show the zero level of each spectrum. The dashed vertical lines 
indicate the limits of each velocity range defined in Table|2] The 
dotted vertical line shows the systemic velocity (-11 km s -1 ). 
The spectral resolution is 3 km s for all lines. 



4. Discussion: physical conditions 

To constrain the physical conditions, we pe rformed radiative 
transfer calculations with the RADEX code (Ivan der Tak et al.l 
2007) based on the large velocity gradient (LVG) approxima- 
tion and assuming a plane-parallel 'slab' geometry for the es- 
cape probability formula. The molecular data were retrieved 
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13-12 
12-11 
2-1 


4.8(0.9) 14.9(1.0) 
6.0(1.0) 19.5(1.2) 
5.1(0.2) 8.8(0.3) 




16.1(0.8) 
17.7(0.9) 
38.9(0.2) 


EHV: -140 to -100 km s 1 ; IHV: - 
SHV: -50 to -20 km s" 1 


-100 to -50 km s _1 



Statistical errors in parenthesis. 




Fig. 3. LVG results produced by the RADEX code for the ob- 
served intensity ratios in the EHV, IHV, and SHV ranges (left, 
middle, and right panels, respectively). The results are shown as 
T vs. n and n x T (thermal pressure) vs. T plots (lower and upper 
panels, respectively). The CO (2-1 )/( 12-1 1) ratio is indicated in 
green and the CO (2-l)/(13-12)ratio in blue. The observed ratio 
is drawn with solid lines, while the line ratio uncertainty (includ- 
ing statistical and 20 % of the calibration errors) with dashed 
lines. N = 2xl0 17 , 5xl0 17 , and 7xl0 17 cirT 3 ; and Av = 40, 50, 
and 30 km s" 1 , for EHV, IHV, and SHV, respectively. 



from the LAMDA data basqj. The c ollisional rate coefficients 
were adopted from Yan g et all (1201 Oh . who calculated the colli- 
sional rates between CO and H2, incorporating energy levels up 
to / = 40 for kinetic temperatures of up to 3000 K. 

For the purposes of our study of the physical conditions as a 
function of velocity, we use the integrated intensities in the EHV, 
IHV, and SHV ranges (Table|2]l. The value that we adopt for the 
linewidth, Av, is directly inferred from the definition of our dif- 
ferent velocity ranges (40, 50, and 30 km s , for EHV, IHV and 
SHV, respectively). Our background radiation field is assumed 
to only be produced by the CMB. Then, our procedure consists 
in running, for each velocity range, a grid of RADEX models to 
compute the integrated intensities within a three-dimensional pa- 
rameter space defined by T (kinetic temperature), n (H2 volume 
density), and N(CO) (CO column density). 

The LVG results were analyzed based on the (2— 1 )/( 1 2— 1 1), 
(2— 1)/(13— 12) line ratios and the absolute integrated intensities 
of those three lines. We point out that owing to the use of these 
ratios the results are biased toward the lower excitations traced 
by the low-/ CO lines, and hence our results should be assumed 
to be lower limits. We use the IR size of the BI (~3"x4") and 
BII (~8"x8") knots to compute the brightness temperature, cor- 



3 http://www.strw.leidenuniv.nl/~moldata/ 
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recting for beam dilution effects. For the SHV component, we 
assumed the size of BII, while for the IHV and EHV compo- 
nents the size of BI was assumed. These assumptions are based 
on the knowledge that the EHV and IHV low-/ CO emission 
peaks around B I, while the SHV emission peaks around BII 
dHatchell et al.lfl999l) . The LVG solutions were found when both 
the ratios and absolute integrated intensities match in the {T, n) 
plane. The N(CO) ranges constrained in this way are 4-9xl0 17 , 
4-6xl0 17 , and l-4xl0 17 cm 2 , for the SHV, IHV, and EHV com- 
ponents, respectively. We tested the influence of the size of the 
emission region on our N(CO) conclusions by varying the as- 
sumed size from a compact emission scenario, where the emis- 
sion spreads over half the BI knot size, to a very extended case in 
which the filling factor is equal to one. Over this range, we found 
that N(CO) varies over up to two orders of magnitude. However, 
this change in N(CO) does not modify significantly the solu- 
tion of the line ratios described below. In Fig. [3] we show the 
RADEX results as plots of T versus (vs.) n and (n x T) vs. T. In 
the latter, (n x T) represents the thermal pressure, and is shown 
because in the LVG analysis the results are not as degenerate as 
for T vs. n (see Fig. 0. The (2-l)/(12-ll) and (2-l)/(13-12) 
ratios yield lower limits to the parameters of T > 100 K, n > 
4.2xl0 4 cm" 3 , and (n x T) > 10 7 K cm" 3 . These lower limits 
indicate that the line ratios actually trace high densities and tem- 
peratures, but the remarkable result of our LVG analysis is that 
the thermal pressure tends to be higher in the IHV range than 
in both the EHV and SHV ranges. In addition, it is quiet a sur- 
prise that both EHV and SHV provide similar constraints, since 
it is usually found that the EHV component (probably related to 
the jet) has higher excitation than the SHV component (which is 
likely to trace the outflow cavity). 

Th e H2O mapping observations at 183 GHz of iLefloch et al.l 
d2011l) . which covered the blue lobe and the central region of 
Cep-E with an angular resolution of ~13", also showed promi- 
nent bullet-like spectral features in the IHV range. The strongest 
emission was found around the central region, while weaker and 
no emission at all was found at BI and BII, respectively. Together 
with additional SiO (8-7) observations, LVG calculations made 
for the gas in the central region yielded the physical conditions 
of T ~ 200 K and n ~10 6 cm" 3 , which are comparable to the 
ones reported here, suggesting that the SiO, H2O, and high-/ 
CO emission arises from gas with similar physical conditions, 
despite the different frequency range and positions probed. A 
correlation between the physical conditions of these molecules 
is an expected result given that the considered formation routes 
of SiO and H2O associate themselves with warm and dense gas. 
However additional observations are required to verify the spa- 
tial correlation between these molecules and to provide more ac- 
curate determinations of the physical conditions. 

Previous ISO observations of CO transitions from (14-13) 
up to (25-26), yielded kinetic temperatures and H2 densities 
in the range 200-1200 K and 4xl0 4 -4xl0 6 cm" 3 , respectively. 
These constraints were obtained by fitting LVG models to the 
total integ rated emission of the s e spectrally unresolved h igh-/ 
CO lines dGiannini et alJl200lb iMoro-Martin et al.ll200ll) . It is 
remarkable that our lower limits are in good agreement with 
the values obtained from ISO, although our results are given 
for three distinct velocity ranges over the same global one. 
However, a more complete LVG analysis is needed to strengthen 
our conclusions with respect to previous ISO studies, and to 
more clearly probe the physical conditions prevailing in the ob- 
served regions. Such an analysis should be based on additional 
observations of both lower- and higher-/ CO transitions with the 
GREAT spectrometer, such as the (11-10) and (16-15) ones, the 



latter would also allow us to make direct comparisons with ISO 
observations. 

Spectrally resolved HerschelfHTFl observations have been 
used to study the emission in different velocity ranges in low- 
mass outflows by measuring the ratios of low- to high-/ CO 
lines, which is an approach that is similar to our analysis. The 
kinetic temperatures generally inferred by those studies f or the 
outflow compo nents are in the range 100-200 K ( Yildiz et al.l 
2010; Bierkeli et al. 201 1). These values are similar to our lower 
limits. As already pointed out, the line ratios of a few high-/ 
CO lines with low-/ CO lines bias the results toward lower ex- 
citation conditions, and could in the end be inappropriate owing 
to the different physical component possibly related to each of 
them, as it has been proven in a few outflows (Ivan Kempen et all 
2010). In this respect, additional observations of more high-/ 
CO lines, as well low- and mid-/ CO transitions, are required to 
test the possibility of the contribution of different physical com- 
ponents to the CO emission in outflows. Eventually, the com- 
bination of a maximum number of CO lines will provide the 
optimal way to shed light on the physical processes responsible 
for the existence o f these components (e.g. iGusdorf et al.ll2012c 
IVisser et al .1120121) . 

5. Conclusions 

The principal conclusions of the present study are: 

1 . Our analysis of SOFIA/GREAT data has demonstrated that 
high-/ CO lines are a good tracer of molecular bullets in 
protostellar outflows. 

2. The bullet at intermediate-to-high velocities has a higher 
level of excitation than the low and extremely high velocity 
gas, at the observed position of the Cep-E outflow. 

3. The still considerable low-velocity emission in high-/ CO 
lines should be taken into account when modeling data. 

4. More and higher-/ CO transitions must be observed to 
break the degeneracy in the LVG solutions. 
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